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Abstract 

We point out a practical way to induce a finite charge Hall conductivity in graphene via an 
off-diagonal strain. Our conclusions are based on a general analysis and supported by numerical 
examples. The interplay between a substrate-induced gap and the strain is discussed. We show 
that the substrate-induced gap itself does not break the symmetry that maps one valley onto the 
other, and hence does not yield a net charge Hall effect. The net charge Hall conductivity is found 
to be sizable and should be observable in the presence of both the gap and the strain. The sign of 
the Hall conductivity can be tuned by varying the strain parameters. The valley Hall conductivity 
is studied as well. The valley and the orbital magnetic moment accumulations on the sample edges 
may be detected in a Hall conductivity experiment. 

PACS numbers: 73.63.-b, 72.80.Vp, 73.22.Pr, 81.05.ue 
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I. INTRODUCTION 



The hallmark of graphene, a single layer of carbon atoms arranged on a honeycomb lattice, 
is its particular band structure with a linear dispersion of the low-energy carriers around 
the two inequivalent (Dirac) points K and K' at the corners of the first Brillouin zone^. 
The K and K' valleys are distinguished by the valley (isospin) degree of freedom which is 
robust against perturbations (unless a momentum on the scale of the graphene's reciprocal 
lattice vector is transferred)-^, which makes it attractive for (valleytronics) applications^"-. 
Hence, it is highly desirable to explore whether transport properties are controllable via the 
valley degree of freedom. As well-established^, the charge Hall conductivity a^y is governed 
by the Berry phase, which for each Dirac cone is finite (±7r for K{K')), but the sum is zero. 
The Berry curvature mediated Hall effect in graphene has been studied in the presence 
of a spin-orbit couplingi^i"— , or for a substrate-induced gap^i^^"— , or in the proximity to 
a superconductor^®. For a monolayer graphene with a uniform interaction gap, e.g. as 
resulting from the interaction with a substrate, one finds^ a finite Berry curvature and a 
finite charge Hall conductivity axy^ ^ in each of the Dirac cones by the interaction gap^^. 
The total charge Hall current vanishes however, for the Berry curvature and the charge 
Hall conductivity for the two Dirac cones are equal with opposite signs. To circumvent this 
problem and obtain a net charge Hall current the authors of Ref.^ proposed the creation 
of a non-equilibrium population in the Dirac cones, which is inherently hard to realize and 
sustain in a real device. 

An energy gap in the electronic structure is advantageous for graphene applications. E.g., 
a zero gap in graphene prevents the graphene field-effect transistor (GFET) device from pos- 
sessing a turn-off state^®. A substrate- induced gap is not restricted to the SiC substrate but 
is also observed for the hexagonal boron nitride substrate^^. A further way to tune the 
energy gap is via the strain engineering, e.g. as brought about by strain super lattices or 
wrinkles^. The strain can be viewed as a symmetry-adapted effective, pseudo magnetic field 
in the sublattice space for the two valleys^"—. A uniaxial strain can be realized experimen- 
tally by bending a flexible substrate^^. And a biaxial strain in graphene can be created by 
shallow depressions^^. Recently, a zero-field quantum Hall effect in graphene was proposed 
via a designed strain profile^®. Usually, strain effect is studied via a minimal coupling which 
is momentum independent's.. Here we go beyond the minimal coupling and include terms 
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that are linear in the momentum by using the theory of invariants^!. The reason is that 
these terms will change the topology of the Dirac cones and which is manifested in real 
physical quantities, for example the Hall conductivity. Furthermore, the effect of strain is 
usually thought to introduce time-reversal-invariant terms in the Hamiltonian. However, 
we show in a time- reversal-breaking environment, for example under a weak magnetic field, 
the strain could couple with graphene via a time-reversal-breaking term which gives rise to 
exotic effects, e.g., under the strain, a net charge Hall current can be driven by an applied 
electric field. This effect may be useful to design strained-graphene-based devices where the 
Hall current can be turned on and off via an in-plane weak magnetic field. 

II. MODEL 

For a general analysis, let us start from a Hamiltonian derived by the theory of invari- 
ants that accounts systematically for all possible strains^. Symmetry operations and their 
compatibility with the underlying lattice structures are encapsulated in this approach. 

The general structure of the two-dimensional (2D) Dirac Hamiltonian is 



The vectors d(k) and d'(k) map the 2D momentum space to a 3D parameter space. (io(k) 
and d'^ik) are structureless parameters in the sublattice space and yield no contributions 
to the Berry curvature and Hall effect (see the appendix for details). The concrete form of 
these vectors depends on the realized physical situations and is given and discussed in full 
details in the appendix, d is parameterizable as d = (i(cos sin6', siny9 sin 6', cos6') (for d', 
we use d', ip' and 6'). The eigen energies are E± = (io(k) ± c?(k) and E'^ = (io(k) ± d'{k.). 
Solving for the eigenf unctions of Eqs. ([1]), we can obtain the Berry connection (i.e. the 
fictitious vector potential as explained in the appendix) and the Berry curvature (i.e. the z 
component of the associated magnetic field in the momentum space) as 



where u = +(— ) for conduction (valence) band. The same applies for Qi^k'z, with 6 and ip 
being replaced by 9' and cp'. The charge current operator is ji = Introducing as the 
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hxi^) = do(k) + cr ■ d(k), for K cone, 
hK'0<.) = c^o(k) + o- ■ d'(k), for K' cone. 



(1) 





Pauli matrix along the z direction to describe the valley degree of freedom we write for the 
valley current operator = \{'^z,ji}, where i = x,y. From the standard Kubo formula^^ 
given in the appendix, we find the charge Hall (cr^y) and the valley Hall {cr^y) conductivities 

as 



^C(V) = ^ K' 

xy xy xy 

2 



(3) 



e 



J2 nMQ,Kz{k) ± J2 n,{k')Q,K'z{k') 

where Aq is the area of the system, n,^(k) is the Fermi function. The Berry curvature can 
be further expressed as 



^uKzik) 



V 



2ci2(k) 



^af^jjaxj /^yd-yf 



(4) 



where is the anti-symmetric tensor, 
a, /3, 7 = 1, 2, 3 and = d^/d^k). 



^^Q^^^ is the pesudo spin current tensors, 
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FIG. 1: (color online) The d vector space which is mapped from a circle in the momentum space 
is shown for the case in absence (a) or presence (b) of strain. The unit of d vector is vpk. The 
other parameters in (b) are |^ = |^ = = = 0.2. 
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III. HAMILTONIAN WITH STRAIN AND ANALYSIS 



A. Analysis for the case without strain 

In absence of the strain and the interaction gap the vector cl(d') has only in-plane com- 
ponents, i.e. di{d'j) = ±vpkx and ^2(1^2) = vpky. Geometrically, a circle in momentum 
space (constant |k|) is mapped onto a circle in the d vector space shown in Fig. [D^a). The d 
and d' are related by a reflection (a reflection through yz plane which is perpendicular 
to the graphene layer see^^). This operation conserves the sublattices but exchanges the 
Dirac cones. Under this operation, [k^, ky) in K cone is transformed to {—k^, ky) in K' cone 
which means the angle of k (i.e. tan0 = ky/k^) is mapped upon Ry^ to an angle ir — (p 
of {—kx, ky) in K' cone. The role of in the d vector space is manifested as a reflection 
through ^2(^3 plane. Therefore, we have ip = (p, (p' = it — (p and 9' = 9 between the phases of 
d and d'. Thus, with increasing in cone (d), the reflected vector in K' cone (d') varies 
in the opposite direction which explains the cancelation of the total charge Hall current (see 
Fig. dl^a)). The above analysis is also valid in the presence of an interaction gap although 
a d{d')^ is introduced and a nontrivial topology in each Dirac cone is induced by it. The 
interaction gap does not break the symmetry between the two Dirac conea^. 



B. Analysis for the case with strain 

To discuss the effect of the strain, let us flrst deflne the strain tensor as 

^ _ 1 _^ ^ _^ /gX 

2 \ dvj dvi dvi dvj J ' 

where u is a polar vector indicating the displacements of atoms. In the theory of invariants^, 
the strain tensor couples with k giving rise to the irreducible tensors K, under D-^h group. We 
only focus on the H^^ which is derived by a multiplication of the irreducible matrices and 
the irreducible tensors that belong both to the Fs representations. Thus, the Hamiltonian is 
expressed in terms of Fi, F2 and Fg, respectively. The corresponding irreducible matrices for 
these representations are 1, cr^, and {ax,<Jy)^- Hence, we infer do = K,^, d^ = /C^, di = K,^'^ 
and d2 = /C^'^ where /C"^''^ indicates irreducible tensors that transform according to the n-th 
irreducible representations of D^h—. 
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From the theory of invariants^! (cf. also the appendix), we obtain the d vector in the K 
and K' cones without including the external electric field 

di = {vfK + V ■ k)r2 - Ax, 
d2 = vpky + (k X v)^ + AyT^, 

ds = f^{vxky + Vykx), (6) 

where denotes the valley degree, vp = VF+h2{uxx + Uyy), Vx = Ks^Uxx-Uyy), Vy = 2be3Uxy, 
and P = 621 /&63 being a dimensionless parameter. As in studies on semiconductors all the 
parameters (6's) are to be determined experimentally or from ab-initio calculations^". The 
effective vector potential whose components are Ax = &6i(^xx — Uyy) and Ay = hQi{2uxy) 
was studied as a minimal coupling in a strain-induced transport^"—. This vector potential 
can also be expressed in terms of the hopping parameters in a tight-binding formalism 
Ax = ^{t?j — ^2) and Ay = \{t2 + t'i — 2ti) where t's are the hopping parameters under 
strain^^. This term is independent of k and can not change the topology of the Dirac 
cones. Please note the third component, i.e. d^, is present under the time-reversal-breaking 
environment. This term violates the d vector from 6 = 7t/2 plane. Thus the d vector 
now describes a map from a 2D momentum space to a 3D vector space. We will study the 
specific consequences of this mapping. Furthermore, the diagonal strain Uxx — Uyy and the 
off-diagonal component Uxy yield new effects, beyond the minimal coupling, and give rise 
to extra terms that depend on k. Therefore, the topology of the cones is changed and a 
non-vanishing Hall conductivity will be the first-order effect from these terms. 

Note, that the term in the above equations stems from the contribution of a smooth 
rippling of the graphene sheet^'^ and results in an isotropic renormalization of the electron 
velocity. The action of the operator is now R~^{Ax, Ay) = {Ax, —Ay), Ry^{vx,Vy) = 
{vx, —Vy). Because of these symmetry relations, in the presence of a strain, a circle in the 
momentum space is mapped onto ellipses in the d vector space, as visualized in Fig. [T](b). 
The two ellipses are reflection-symmetric at did^ plane in the d vector space. This plane is 
perpendicular to the one in the case of zero strain. This is the consequence of the symmetry 
property for the strain tensor Uij under Ry^- Uij is not a polar vector (say k) and transforms 
like the symmetrized {ki, kj}. However, the symmetry plane of did^ links d(k) and d'(— k) 
and vice versa (rather then d(k) and d'(k)) (see Fig. [It^b)). We note d^lk) = —d'^^—k). 
Therefore, we can rearrange the summation in deriving a^y^'' over k' in K' cone by the 
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restriction k' = — k in Eq. ([3]). So the phases of d' are 0' = 27r — and 6' = tt — 6. According 
to Eq. ([2]) we obtain the remarkable result 

fluKz = ^uK'z- (7) 

The two Dirac cones have the same conductivity contributions in contrast to the opposite 
contributions in zero strain case. The opposite properties of the Berry curvature for the 
conduction and the valence bands are maintained since the strain does not break the particle- 
hole symmetry. Eq. ([7]) evidences a net charge Hall conductivity in our case. 

IV. NUMERICAL ILLUSTRATIONS 

To demonstrate the symmetry analysis, in Fig. [2] we show the results of the numerical 
calculations for the total Berry curvature in (a), (b) and (c). The observations from (a) and 
(b) are: 

i) the prominent peaks (or dips) indicate the positions of the Dirac points in momentum 
space shifted by the effective magnetic field induced by the strain. 

a) A singularity of these peaks (or dips) does not exist generally since a strain- induced gap 
is present. 

in) The Berry curvature apart from the Dirac points is generally finite. 

iv) The two peaks in Berry's curvature are just at the antipodal points with respect to the 
zero k (inversion in the momentum space). 

v) The most important fact is that the Berry curvatures at these two points possess the 
same sign. These observations are exactly what one expects from a symmetry analysis. 
Furthermore, we find the Berry curvature is exactly zero if the off-diagonal component Uxy 
is zero. For Uxy 7^ 0, as shown in Fig. Mjo), the Berry curvature is non vanishing even when 
Uxx = Uyy (which means Ax = Vx = Q)- Note, the discussed symmetry is still preserved with 
the result of a net Hall conductivity. The difference to Fig. [2]^a) and (b) is that the Hall 
conductivity is present with a reversed sign. The energies (with respect to (io('^o)) for the 
conduction and valence bands, and the total Berry curvature calculated along the dash line 
in (a), are shown in Fig. WA) and (e), respectively. As evident, the two peaks in the Berry 
curvature are the same and appear at the positions where the smallest band gap exists. The 
inset in (e) shows a zoom-in of the opened strain-induced gap. 
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FIG. 2: (color online) The contour plot of the total Berry curvature r^tot for the conduction band 
is calculated numerically in (a), (b) and (c). The energy profiles in (d) and (f) and the Q.tot in (e) 
and (g) are calculated along the dashed lines in (a) and (c) respectively. The inset of (d) and (e) 
is a zoom-in plot around K' cone. The units are = 1 nm~^, vp^ vpfikQ for the wave vector, the 
velocity, and energy respectively (note the unit of A is energy). /3 = 0.5 is kept for all graphs. A 
energy gap A = 0.28 eV— is present in (c), (f) and (g) and zero elsewhere. The other parameters 



Ay 



vp vpko vp vpko 



0.3 in (a), (c)-(g); 1^ = ^ = 0, |^ 



Ay 



Vp Vpko 



0.3 in (b). 



In Fig. |2]^c) shows the calculated the total Berry curvature for the conduction band in 
the presence of a strain and an interaction gap^. The interaction gap is uniformly generated 
for both Dirac cones- and thus adds to d-^ and d'^. The interaction gap breaks the symmetry 
shown in Fig. Wljo). If (is = = A/2 only, the odd-symmetry of the Berry curvature for 
the two cones is preserved. If there is only a strain effect in d^ and d'^, the symmetry is 
changed to even. If the interaction gap and the strain are both present, they compete with 
each other. Generally, the Berry curvature does not possess an explicit symmetry which 
can be seen from Figs. [2]^f) and (g) along the dash line in (c). The gaps at the shifted K 
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and K' cones composed of the interaction gap and the strain-induced gap are generally not 
symmetric now. 
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FIG. 3: (color online) The dependencies of a'^y on /i in (a) and on and Vy in (b) and (c). (d) 



shows the variation of a^y with Vx and Vy. In (a), Vx = Vy, a = 1 and /3 = 0.5. The parameters in 
(b), (c) and (d) are /3 = 0.5, a = 1.2 and ^ = 1.0. The unit of the conductivity is /h. The other 
units are the same as in Fig. [2j 



Introducing the relation Ax{y) = OiVx{y) we calculated the total charge Hall conductivity 
with varying the Fermi level (Fig. [3] (a)). For graphene without the interaction gap there 
is no sign change in a'^y. In the presence of the interaction gap and the strain, there is a sign 
change since the interaction gap is modified to two non-symmetric gaps in the two cones 
(see Fig. |2](c) and (f)). With increasing /i, one cone dominates first. The other cone follows 
up with a further increase of leading thus to a sign change. The values at which the sign 
change occurs indicate the asymmetry of the gaps at K and K' . The saturation value of the 
a'^y is increased with a larger strain. To obtain a deeper insight into how the aSj^^ varies 
with the strain parameters, the contour plots are shown in Fig. [3t^b)-(d). In (b) and (c), a^y 
is shown for A = and A = 0.28 eV, respectively. The regions of a finite Hall conductivity 
are split into islands by a zero regions. The dash lines in (b) and (c) indicate the relation 
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Vy = \/3vx from which d{d') = can be obtained. When this happens there is no structure 
in the system. In the upper part of this hne, the islands show a positive sign of the charge 
Hall conductivity. In contrast, in the lower part of the line, the charge Hall conductivity 
is negative. In regions (HI) and (IV) in Figs. [3]^b) and (c), the strain effect is prominent 
leading to a large charge Hall conductivity. The areas of the regions (HI) and (IV) shrink 
accordingly in Fig. ^c) with nonzero interaction gap. An interesting observation is that 
a transition may occur when crossing from the lower part to the upper part. Noting that 
the strain is coupled with the graphene in a symmetry-governed way, the symmetry at the 
transition point is maintained and the topology is changed. This could be an example of a 
Lifshitz transition^^ induced by the strain in graphene. In the regions (I) in Figs. [3]^b) and 
(c), a^y is small since the interaction gap dominates the strain- induced gap in (c) and a small 
strain-induced gap in (b). The regions (II) in (b) and (c) evidence that the induced gap is 
large leading to a zero for each cone. The sign change of is caused by the reversed 



sign of the cross products of the two pseudo spin currents, i.e. ix{y) = {jix{y), j2x{y), jsxiy)) in 
3D d space. Fig. |3]^d) shows the valley Hall conductivity . In contrast to the corresponding 
regions in (b) and (c), the a^y is large in region (I) and possesses the same sign for the (HI) 
and (IV) regions. In regions (II), the valley Hall conductivity is quite small as well. 
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FIG. 4: (color online) Schematics of the charge and valley Hall effect under an external electric 
field for a zero strain and nonzero interaction gap in (a), nonzero strain and zero interaction gap 
in (b), and nonzero strain and nonzero interaction gap in (c). The arrowed curves denote the 
electron flows. The arrowed circles indicate the orbital magnetic momentum in K and K' cones. 
The thicker (thinner) curve means larger (weaker) electron current. 
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V. POSSIBLE EXPERIMENTAL SETUP 



Experimentally, the valley degree could be characterized by the valley-carried orbital 
magnetic moment^i^ii^. In our study, the orbital magnetic moment in the presence of the 
strain is rriuz^ = — z/2^^fi,^x(x')z(k), where z/ = ± for the conduction and the valence 
band. This formula reduces to that in absence of the strain given in Ref.^ except for a minus 
sign because of a sign change in the definition of the Berry curvature. As is known, the 
orbital magnetic moment for a Bloch wave packet stems from its self rotation and can be 
tested in an external magnetic field. 

For a zero strain, no charge Hall current is present, as illustrated in Fig. lU^a). The 
electron flows from K and K' cones carry different orbital magnetic moments, which gives 
rise to a finite valley Hall conductivity. The opposite orbital magnetic moment accumulates 
on the opposite edges of the Hall bar and may be tested in experiments, e.g. by giant 
magneto resistent sensors. For nonzero strain and zero interaction gap, the net charge Hall 
current is nonzero since the electrons fiow into the same direction with the same orbital 
magnetic moments in the two cones (see Fig. 111(b)). The interesting point is that the net 
magnetic moment is non vanishing as a consequence. For the nonzero strain and nonzero 
interaction gap, the electron fiows may proceed in the same directions; or they may run 
in opposite directions with different magnitudes of current density and orbital magnetic 
moment. Which case is realized in determined by the relative ratio of the strengthes of 
the strain and the interaction gap (the latter case is schematically shown in Fig. Hl^c)). 
The effect could be observed by applying a weak magnetic field in the plane (to minimize 
the orbital motion and the Zeeman effect could be negligible) to the graphene device with 
strain or by depositing the graphene device on a (magnetic) substrate which breaks the 
time-reversal symmetry. 

In summary, we inspected the effects of strain in graphene within a model that goes 
beyond the minimal coupling by utilizing the theory of invariants. We show how the strain 
tunes the topology of the two Dirac cones, the associated Berry curvatures, the orbital 
magnetic moments and the conductivities. Unlike the other studies, we also consider the 
effect of strain under a time-reversal-breaking environment and point out that a net charge 
Hall effect could be induced by an off-diagonal strain. This is a case which falls in the 
category of the anomalous Hall effects without spin orbit coupling and could be tested by 
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applying a weak magnetic field in the plane of the graphene device and may be useful in 
new potential applications in strained-based graphene devices. 

Acknowledgements: We thank R. Winkler and V. K. Dugaev for useful discussions. 
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Appendix A: The Hamiltonian in presence of a strain derived from the theory of 
invariants 




In Fig. El^a), the graphene lattice structure is shown, oi and a2 are the primitive vectors 
of a unit cell. Fig. [5|^b) shows the corresponding first Brillouin zone. The crystal point group 
at the r point is D^^ and at the Dirac points, i.e. K and K' , is -D3/1. We do not present 
the detailed group representations here explicitly, for those have been discussed in depth in 
Refs.— "2^. The key point is that -D3/1 is a subgroup of D^h. Restricting the considerations 
to the low energy regime around the two Dirac points, the constructed Hamiltonian has to 
be invariant under the symmetry operations in D^^. Therefore, the symmetry operations 
belonging to -Dg/i, but outside D-^h will not preserve the invariance of the Hamiltonians at the 
points K and K' . These two Hamiltonians at K and K' are interrelated via the symmetry 
operations in D^h group which map K onto K' and vice versa. Particularly, we consider one 
such operation Ry which is a reflection plane through yz (and kykz in the momentum space) 
perpendicular to the graphene layer. Under this operation, the sublattices are reserved but 
K is mapped onto K' . We then have 

/iK'(k) = V{Ry)hK{Ry'\<i)V-\Ry), (Al) 

where V{Ry) is a representation matrix for the operation Ry. 
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Let us define the strain tensor as 

1 f dui duj duzduz\ 

2 \ dvj dvi dvi dvj J ' 

wliere u is a polar vector indicating the displacements of atoms. In the theory of invariants^, 
the strain tensor couples with k to give rise to the irreducible tensors /Cs under D^h group. 
We only focus on the H^^ which is derived by a multiplication of the irreducible matrices and 
the irreducible tensors which both belong to the Fs representations. Thus, the Hamiltonian 
is expressed in terms of Fi, F2 and Fg, respectively. The corresponding irreducible matrices 
for these representations are 1, o"^, and (cr^, UyY'^ ■ These four matrices form a complete space 
in which any 2x2 quantity can be expanded in terms of them. The expansion coefficients 
are just the irreducible tensors /Cs. In our concrete study about the strain effect, the d 
vectors in Eqs. (1) and (6) in the main text are related to the irreducible tensors, = K}, 
d^ = K? , di = /C^'^ and ^2 = /C^'^ where /C'^ indicates irreducible tensors that transform 
according to the K-th irreducible representation of D^^h- These irreducible tensors except the 
one for Fi span the d vector space in K cone. Therefore, we can construct a Hamiltonian 
in the K cone which is invariant upon D^h group 

/i;^(k) =rfo(k) + o--d(k). (A3) 

Explicitly, the d vector in the K cone reads 

di = vfK + V ■ k - a, 
d2 = vpky + (k X v)^ + Ay, 

ds = (3{vxky + Vyk^), (A4) 

where vp = vp + &62(mxx + Uyy), = bQ3{uxx - Uyy), Vy = 2bQ3U^y, A:r = &6i(mxx " Uyy) and 
^y = bei{2uxy), and (3 = 621/^63 being a dimensionless parameter. The parameters bs are 
material-specific parameters which need to be determined by experiments. The (io(k) will 
not make a contribution to the Hall effect and Berry's curvature as shown below. However, 
we write it exphcitly here for a reference, i.e. do(k) = 611 1 + bi2{uxx + Uyy)l + &i3[(%y — 
Uxx)kx + 2uxyky\l. For zero strain, the first term is related to the chemical potential. It will 
be zero for half-filling. 

From Eq. (]A1|) and note Ry^{Ax, Ay) = {Ax, —Ay), Ry^{vx,Vy) = {vx,—Vy), we obtain 



14 



the Hamiltonian in the K' cone 



hK\\^) = rf'o(k) + cr.d'(k), 
d'l = -vfK - V • k - 
d'2 = vpky + (k X v)^ — Ay, 

4 = ds. (A5) 

Introducing r^, the third Pauh matrix, to indicate the valley degree of freedom, its eigenvalue 
is ±1 for K{K') cone. Therefore, Eqs. f lASP and f lASP can be compacted into the Hamiltonian 
given by Eq. (6) in the main text. 



Appendix B: Berry's curvature 

For the Hamiltonian Hk, we parameterize d as d = (i(cos(y9 sin ^, sin sin^, cos ^) (for d', 
we use d', ip' and 9'). The eigen energies are E± = do(k)±d(k) (for h^', E'^ = d'Q{k)±d'(k)). 
The eigenvectors for hx are 

T 



6 6 
\^+k) = [ cos-, sin -e^'^ 

..^ ^ 



\^-k)= (^sin-,-cos-e^'^J , (Bl) 

where ± stand for the bands with the energies E±. Similar equations for the K' cone can 
be derived by substituting 0',ip' into the corresponding places in Eq. (1B1|) . The Berry 
connection (i.e. a k-space vector potential) in the band u is given 

a,,(k) = -t{^,r\V^\^,r), (B2) 

where r being the valley index. The corresponding Berry curvature (a k-space magnetic 
field) is 

n,r(k) = Vk X = n,rzi'k)e,. (B3) 

For graphene, there is only the component in z direction. Therefore, we can obtain the 
Berry curvature by using the eigenvectors as'^- 

sin 0^ ( 80 r dipr d0rdipr\ .-R.x 
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1. Zero strain and zero interaction gap 

We have hK(K') = ^o'xdi + ayd2, di = vpkx = vphcoscf), d2 = vpky = vphsmcj), and 
^3 = 0. Thus, we infer = (p and = ^ for K cone. It is straightforward to show that 

1 , . 

= ^60. (B5) 

This is a vector potential of a monopole in the momentum space (note, there is sign difference 
in Eq. (IB 5 1) ("-" instead of "+") to that in Ref.— because of a sign difference in the 
definition). The Berry curvature is given as = (l/2)<^^(k)- The Berry curvature 

changes its sign from one band or valley to the other. It is now evident that there is a 
singularity at the Dirac points in the momentum space. 



2. Zero strain and nonzero interaction gap 

We have di = vph^, d2 = vpky, and ^3 = -y. Thus, we have (f = (j) and 6 = tan~^ for 
K cone. Straightforward computations deliver for the Berry connection a+x = — ^—ff^e^ 
and for the Berry curvature 

3AaH^ 
^""^ ~ 2(A2 + 3/tW)3/2' 

where vp = ^o,t is used. These results match exactly those derived in Ref.- with a sign 
change from "+" to "-" because of the difference in the definition (the definition in our 
paper is chosen so that the results are consistent (including the signs) with the results of 
the Kubo formula). 



3. Derivation of Equation (4) 

For a general 2D Hamiltonian for a Dirac fermion, like Eq. ( 1A3I) . we can calculate the 
charge Hall conductivity driven by an external electric field. The current density operator for 
the K cone is J/^(k) = = joi + jaiO"", where a = 1, 2, 3, i = x, and jai = ^^^r^- From 
tan^ = ■\/d\ + d^/d^, we get ^ = cos*/? cos 6'/(i, ^ = sin cos 6'/(i, and ^ = —sinO/d. 
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Thus, we infer 



d9 09 89 89 



8kx 8di 8d2 8d; 



3 

COS if COS 9 . sin Lp cos 9 , sin 9 
jix 2 ^ -^^x 2 -^sx —J— 

„ I \ ,^ \<i\{3\xdz - j-sxdi) + d2{j2xd^ - hxd2)] 

^ 'dx(j, xd)]3, (B7) 



d'^/dl + dl 

where j^(j^) = {jix{y), 32x{y), hx{y))- Similarly, we find 

89 1 



8ky d^^^T^, 



8K dl + dl^^'^'' 



d X (jy X d)] 



3 ' 



8^p 



"h X . (B8) 



8ky d\ + d^ 

Substituting Eqs. ( ]B7p and (IBSP into Eq. ( 1B4p . we get the Berry curvature 

^ irf3(rffTrfi) ^ ^ ^^^^ ^•'^ ^ ~ ^ ^"^'^ ^ ^^^^ ^"^^ ^ ^^^^ 

= 2^[^°/37i"a:j/3yC?7]- (B9) 

This is just Eq. (4) in the main text. A similar equation for the K' cone can be derived 
by changing the corresponding d and ]x(y) vectors. This formula for the Berry curvature is 
used in the numerical calculations to produce the results discussed in the main text. 



Appendix C: Anomalous charge Hall conductivity 

1. Anomalous Charge Hall conductivity for one cone 

We first calculate the charge Hall conductivity for the K cone. The Green's function 
reads 

GK%iUn) = [iuJn " /ix(k)]"\ (CI) 
where Un are the Matsubara frequencies. This Green's function can be further evaluated as 
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where t = ± and P± = (1 ± daCr°')/2. We use the formuL 



kn 



^>r- Tr [jf (k)P.(k) Jf (k)P,(k)] 

/3^o 4- [^(w„ + i^m) - (k)] - Pf (k)] ' 



kn s,t=it 



(C3) 



here /3 = {ksT) ^ indicates the temperature. In terms of Qf^ we can calculate the charge 
Hall conductivity 

< = ^ lim-Q,^(c. + z5) 



i 



e2 , Tr [Jf (k)P,(k)Jf (k)P,(k)] 

2- — ^^^^ ^""'^^^ ~ ""'^^^^^ 



[P,(k)-P.(k)]^ 



St k 



where P^ = P±, and j/^'^ = Jf-. Thus 

K(k)-n+(k)] 



0" 



K 

xy 



E 



4d2(k) 



X 



Tr 



i 



K(k)-n+(k)] 



Tr 



(jox + jaxCr")(jos, + 3i,y(y^)dy(y'' 



- {jox + jaxC^°)dl3Cr'^{joy + 3t,y(^^ 

\n_ — n 



1 



2rf2(k) 



h 4Ao ^ d^iV) 

k 



-Tr (jox+Jaxa")j^j,rf.K,a^ 
Tr (jox + 



^ An 



(C5) 



!/,k 



where fi^x^lk) = ^£ai3-yjaxji3yd"y IS the Berry curvature. Eq. fICSp is just the charge Hall 
conductivity for the K cone, as given by Eq. (3) in the main text (cf. also Ref.— ). 
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2. The total charge Hall conductivity and the valley Hall conductivity 



Accounting for the valley degrees of freedom, we may utilize Pauli matrices rs to char- 
acterize the valley space. do(k) and rfo(k) can be put into a matrix in the valley space 



rfo(k) 



Do 



The d vector is now D, and we have 



d'M 



H{k) = Do{k) + a--D{k), 



where 



D2 = Vpky + (k X V)^ + AyTz, 

D3 = I3{vxky + Vykx). 
Therefore, the Hamiltonian in the valley space is now 



The charge current density operator is 



/ \ 

hK{k) 



hK'ik) 



dH{k) 
dki 

( \ 




J, 



K> 



(C6) 



(C7) 



(C8) 



(C9) 



(CIO) 
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The valley current density operator is defined as 



( 



dH{k) \ 
\ 







V 



-J, 



K' 



(Cll) 



The Green's function is 



G(k, iun) = 



Prom the Kubo formula it follows 



Q'iPiivm) = ^E^' [jf(^)(k)G(k,iK + ii.J)jf (k)G(k,ia;„)J , 

^ k,n 



GK{k,iUn) 

GK'{k,i(^n) 



(C12) 



(C13) 



where the prime in the Tr' means that we also need to calculate the trace over the valley 
space in addition to the sublattice space. Let us evaluate the trace over the valley space 
first. This yields 



IV [jf (k)G;,,(k, i{Un + iVm))jf{\^)GK'{\^, iu^)] } . (C14) 



Therefore, we obtain 



■xy 



= ^xy^^xy 



(C15) 



The "+" indicates the charge Hall conductivity with the contributions from the two cones, 
and the "-" denotes the valley Hall conductivity. 
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